Effect of magnetic fluctuations on the normal state properties of Sr2Ru04 
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We investigate the normal state transport properties of Sr2Ru04 and we show that a consistent 
explanation of the experimental results can be obtained assuming that the system is near a quantum 
phase transition. Within the framework of a self-consistent spin fluctuation theory, we calculate the 
temperature variation of some relevant physical quantities and we discuss a possible microscopic 
origin of the quantum phase transition. 



Since the discovery of superconductivity in the lay- 
ered perovskite oxide Sr2Ru04 [Q much attention has 
been paid to the electronic, magnetic and superconduct- 
ing properties of this compound. In its stochiometric 
composition Sr2Ru04 is metallic, shows a Fermi liquid 
behavior below a certain temperature, undergoes to a 
superconducting transition at the rather low transition 
temperature of ~1K. 

Recent experiments lead to speculate that the electrons 
in the Cooper pairs would have aligned spins indicating 
a triplet superconductivity, and an odd orbital wavefunc- 
tion as for example p-wave in contrast with the 

spin-singlet pairing of conventional superconductors as 
well as of the high-Tc copper oxides. 

Its electronic properties are mainly determined by 
three bands crossing the Fermi level produced by anti- 
bonding Ru 4de (xy,yz,xz) and O 2p7r orbitals 

It is found that the electrons belonging to the xy band 
play an important role determining the normal state 
properties. From ^^O NMR in ^^O isotope-enriched mea- 
surements 0, js], spin correlations in Ru d^y and dxz.,yz 
orbitals are probed separately. They provide the evi- 
dence that only the spin susceptibility Xxy from the elec- 
trons in the dxy orbital shows significant temperature de- 
pendence, increases monotonically with decreasing 
temperature down to about 40K, turns over, and then 
tends to level below TpL =25K in the canonical Fermi 
liquid state, implying that the spin correlations are pre- 
dominantly ferromagnetic in origin, and orbital depen- 
dent. We point out that, below TpL the out-of-plane and 
in-plane resistivity have a quadratic temperature depen- 
dence indicating that a crossover to a Fermi-liquid state 
takes place below 25K prior to the superconducting tran- 
sition. On the other hand, above Tfl the in-plane resis- 
tivity show a superlinear behaviour which may be related 
to scattering due to spin fluctuations. These experimen- 
tal results indicate a clear connection to the crossover in 
the magnetic transport measurements. 

The demonstration that spin correlations in the Ru02 
planes are predominantly ferromagnetic comes also from 
measurements of nuclear spin-lattice relaxation rate. |^] 
The hyperfine form factor for the ^^^l/TiT at the Ru 
site is costant in the reciprocal space, while is propor- 
tional to cos{qx,y/2)'^ for ^''l/TiT at the planar oxygen 



O site. This means that the ^"^l/TiT can probe both fer- 
romagnetic and antiferromagnetic spin fluctuations. On 
the other hand, since the form factor for the planar O 
site is zero at the staggered wave vector q — {7r,7r}, 
antiferromagnetic spin fluctuations do not contribute to 
Hence, one can probe whether the spin cor- 
relations are ferromagnetic or antiferromagnetic by com- 
paring the i^H/TiT and "l/TiT. It turns out |§ that 
both ^"^l/TiT and increase monotonically down 

to TpL, and almost saturate in a Korringa-like behaviour. 

Besides, inelastic neutron scattering measurements 
|]lO| in the normal state reveal the existence of incom- 
mensurate magnetic spin fluctuations located at q = 
(±0.67r/a, ±0.67r/a, 0) due to the pronounced nesting 
properties of the almost one-dimensional d^z.yz bands. 
Infact, the ID sheets can be schematically described by 
parallel planes separated by g = ±27r/3a, running both 
in the x and in the y directions which give rise to dy- 
namical nesting effects at the wave vectors k = (q,ky), 
k = {kx,q), and in particular at q = (g, q). [|ll| 

Combining the results of NMR, |Q nuclear spin-lattice 
relaxation rate, and inelastic neutron scattering mea- 
surements [0, it is possible to draw the following phys- 
ical picture for the Sr2Ru04. There is a strong en- 
hancement of ferromagnetic spin correlations above Tp^ 
mainly due to ferromagnetic correlations between the 
electrons in the dxy band as revealed from Knight shifts 
experiments. These contributions are in competition 
with incommensurate spin fluctuations due to the nest- 
ing properties of the almost ID dxz,yz bands. Hence, the 
main features of the magnetic response turn out to be 
decoupled for the electrons in the dxy and in the dxz,yz 
bands, respectively. 

The aim of this Letter is to propose an explanation 
of the non Fermi liquid behaviour observed in the Xxy 
and in the 1/TiT experiments of the Sr2Ru04 assuming 
that this compound is close to a 2D ferromagnetic quan- 
tum critical point. Then, we will discuss how it might 
be possible to drive the system from a magnetic to a 
non-magnetic state considering the interplay between the 
lowering of the crystal symmetry and the magnetic cor- 
relations. Let consider the effect of ferromagnetic fluctu- 
ations in the correlated dxy band considering the system 
as a 2D ferromagnet with vanishing critical temperature. 
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FIG. 1. Magnetic susceptibility for the dxy orbital. The 
circles stand for the experimental results from Q and the 
solid line is the theoretical curve. 



The starting point in our case is represented by the corre- 
lated dynamics in the Ru02 planes of the d^^y electrons, 
which we consider responsible for ferromagnetic fluctu- 
ations as suggested from NMR and nuclear spin-lattice 
relaxation rate. 

In this case as in many others metallic systems at low 
temperatures it has been observed the signature 

for a non-Fermi-liquid behavior near a quantum phase 
transition. The proximity to a quantum critical point 
is usually marked by enhancement in the effective mass, 
and in spin and charge response at low temperatures. 
This in turn introduces a low energy scale which marks 
a crossover from quantum to a classical behavior in the 
temperature dependence of various physical properties. 

Now, to calculate some relevant physical properties, we 
consider the self-consistent spin fluctuation theory devel- 
oped by Ramakrishnan et al and many others 
for the dxy band. 

The main motivation for the use of spin-fluctuation 
theory is the largeness of the Stoner enhancement factor 
l/a{0) a,t T — 0. In such a case, a highly paramagnetic 
system at low temperature can be considered to be in the 
vicinity of a magnetic transition. The variation of rele- 
vant physical quantities is therefore governed by trans- 
verse and longitudinal spin fluctuations. Even though the 
order parameter is zero above the transition, the effect 
of fluctuations is observable well above the transition. 

From a general point of view, the susceptibility can 
be obtained from the equation of state for the external 
magnetic field H and the magnetization M 



H 
M 



a{T) + PM 



f2 



(1) 



where the coeflicient a and /3 took different expressions 
depending on the approximation. In the mean fleld the- 
ory of itinerant ferromagnets aupiT) = 1 — Up{eF) and 
f3 is constant, pi^p) being the density of states at the 
Fermi level. As the temperature dependence is weak, 
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FIG. 2. Nuclear spin-lattice relaxation rate at Ru site vs 
temperature: the solid line stands for the theoretical result 
and the circles are the experimental data taken from Q. 

of the order of {T/Tp)'^, one cannot yield a Curie- Weiss 
form of the spin-susceptibility. Here Tp- is the Fermi tem- 
perature. This issue is overcome in the spin fluctuation 



theory where a{T) is given by |13 14[ |: 

a{T) = a(0) + u{2D^ + 3D^) 



(2) 



here a(0) is the susceptibility enhancement factor at 
T — 0, which includes the mean field part of Eq. |l] and 
the zero-temperature part of the fluctuation self-energy. 
Still, and are the transverse and longitudinal 
spin fluctuation amplitudes obtained by a diagrammatic 
expansion of the effective Hamiltonian, after integrating 
over the electronic degrees of freedom. The factor u is a 
dimensionless short range fluctuation coupling constant 
|l3| . In this framework, the susceptibility for the ferro- 
magnetic case has the following form: 



a(T) + (5q2 



(3) 



where the parameters 6 and 7 depend on the diagram- 
matic expansion. 

It is found that there are two different regimes which 
for a 2D fcrromagnet read as follows: when r < a{Q), one 
has an enhanced Pauli susceptibility behaviour, while at 
higher temperatures, i.e. for a(0) < r < 1, the suscep- 
tibility follows the classical Curie Weiss behavior. Here, 
T=T/Ti?. Therefore, the quantity a{0)Tp plays the role 
of the low energy scale and naturally comes out from the 
fluctuation theory. The comparison between the exper- 
imental data from Imai et al [Q and the xy spin sus- 
ceptibility calculated within the above assumptions, is 
reported in Figj^. As we can see, there is a good agree- 
ment between the experimental results and the theoreti- 
cal prediction. It is worth stressing that using the exper- 
imental data for the spin magnetic susceptibility and the 
Fermi temperature, calculated using the Fermi wavevec- 
tor and the effective mass reported by Mackenzie et al 
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FIG. 3. In plane resistivity vs temperature: the circles are 
the experimental results from jisj while the solid line is the 
theoretical curve. 

we obtain for a{0)Tp the value ^70K which gives 
the right order of magnitude for the temperature T fl of 
the crossover into a Fermi liquid state. 

Within the same approach it is easy to obtain the 
temperature dependence of the nuclear spin-lattice re- 
laxation rate. 

One gets {TiT)-^ - a[Tf/'^ for a normal Fermi- 
liquid a{T) is constant and gives the usual Korringa law, 
while for a 2D ferromagnet varies as T~^. 

Assuming the previous temperature dependence for 
(TiT)-i for the xy electrons and a usual Korringa law 
for the other two bands, we obtain a fair fit, shown in 
Fig.^, to the experimental results of Imai et al Q for 
temperatures larger than T^fl- 

It is worth pointing out that though the contributions 
of the dxz.yz contain the effect of dynamical nesting, in 
the present calculation they have been added at a mean 
field level so that there is a weak temperature dependence 
of the order of ~ (T/Tp)^ in the range examined. If we 
consider the correction due to spin fluctuations instead 
that at q = around the nesting vector q = (g, q) we 
would expect a contribution to (TiT)~^ with a temper- 
ature dependence ^ T^^ . A more detailed analysis con- 
sidering the multi-orbital system within the spin fluctu- 
ation theory is in progress. Furthermore, we expect that 
the transport properties are affected by spin-fluctuation 
above Tfl and that this effect should manifest mainly 
in the in-plane quantities. Therefore, we have calculated 
the in-plane resistivity using the Boltzmann equation for 
a multi-band anisotropic electronic system assuming for 
the relaxation rates, two Fermi liquid {tJ^^ ~ T^) behav- 
iors and a two-dimensional ferromagnetic decaying time 
^ T as due to the scattering of the spin fluctuations 
with the electrons in the dxy band [p^ , p^ . 

It is worth pointing out that the relaxation rates Tf 
have been determined in a way to reproduce the in-plane 
Hall coefficient and the magnetoresistance, respectively 
[p7t . The fit to the experimental data reported in FigJ|, 



where the experimental data are taken from [|18|, indi- 
cates that the assumption of linear temperature depen- 
dence in the decaying relaxation rate of the xy electron- 
dynamic is essentially correct. 

Let consider now the origin of the ferromagnetic spin 
fluctuations looking at the problem from a more general 
point of view. 

One argument on why the ferromagnetic spin fluctua- 
tions play an important role in determining the transport 
properties of Sr2Ru04 is the following. The Ruddlesen- 
Popper series Sr„+iRu„03„+i contains multilayer com- 
pounds with n as the number of Ru02-planes per unit 
cell. In particular, there is a series of ferromagnetic 
compound related to Sr2Ru04: the infinite- layer (3D) 
SrRuOs is ferromagnetic with Tc ~165K for n=3 

one finds T^ ~148K ||§ and for n=2 Te ~104K 
although the ferromagnetism in this last case is still con- 
troversial. These values indicate the tendency that with 
decreasing the layer number n, Tc is reduced and finally 
vanishes , since for n = 1 the system is paramagnetic 
at all temperatures before becoming superconductor be- 
low - IK. 

Since n is varying in a discrete way, it is meaningless 
to consider it as the direct parameter whose continuous 
change controls the transition from a paramagnetic to 
a ferromagnetic phase. On the other hand, it is possi- 
ble to image stacks of multilayer Sr„+iRu„03„+i com- 
pounds with different n so that the final system has an 
average number of planes which may be varied continu- 
ously. In this framework, the average number of Ru02 
planes can represent the parameter controlling a quan- 
tum phase transition. This kind of ideal system can in 
principle be realized by performing a molecular beam epi- 
taxial growth. 

An important point on the changes induced by n is 
that when we move from the cubic SrRuOa compound to 
the tetragonal Sr2Ru04, the modification of the magnetic 
ground state is accompanied by the lowering of the lattice 
symmetry. 

We know other examples within the metal transition 
oxides whose magnetic properties are strictly related to 
structural changes. |p3| These cases are mainly insula- 
tors systems where the Jahn- Teller effect cause structural 
phase transitions, and gives rise to orbital ordering and 
new magnetic structure. 

In our situation, the removal of the degeneracy due to 
crystal field effects without Jahn- Teller distortions pro- 
duces a subtle change in the magnetic properties. Due 
to the peculiar filling of the 4d level in the SrRuO sys- 
tems, we will show that the ground state switches from 
a paramagnetic to a ferromagnetic configuration depend- 
ing on the way one lowers the lattice symmetry, i.e. by 
elongating or compressing the RuOg octaedra along the 
c-axis. 

This effect might be experimentally fulfilled by apply- 
ing an external uniaxial pressure and/or by means of sub- 
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stitution of the Sr+^ with a larger ion with the same va- 
lence. In this way, it might be possible to induce changes 
in the local crystal field driving the system from a mag- 
netic state to a non-magnetic one or vice versa. 

To have a quantitative insight on these considerations, 
we study the effect of the removal of the cubic symmetry 
for a cP system (two holes in the t2g orbitals). 

The Hamiltonian is given by: 

i.a 

ia i,a<l3 

where the Hamiltonian contains intrasite terms, and the 
diagonal hopping between electrons in the different t2g 
orbitals. 

The operator da,ia destroys an electron with spin u in the 
orbital a on site i, n.i^aa is the electron density with spin 
a on the orbital a, and Si^Q, is the spin of the electron in 
the a orbital on site i, respectively. Here, Ci^a denotes the 
on-site energy of the a i2g orbital on the site i, U and Jh 
stands for the Coulomb and Hund's exchange interaction, 
respectively, and t denotes the hopping amplitude. 

The cubic-tetragonal symmetry crossover is controlled 
by means of the zero energy splitting A = {f.[dxy) ~ 
e{djz)) between the d^y and dxz,yz orbitals. 

In particular, the variations of A describe the interplay 
between ionic and covalent effects induced by the oxygens 
surrounding the Ru ion and thus simulate the effect of 
compressing and/or elongating the RuOg octaedra along 
the c-axis. 

The phase diagram is obtained considering two effec- 
tive ruthenium atoms. The results of this numerical sim- 
ulation are reported in Fig.^ with Jh /U vs A/t. 

For the cubic case, i.e. when A=0, the ground state 
is ferrimagnetic (total spin is equal to one) in a very 
small region, above that it becomes ferromagnetic pass- 
ing through a spin zero (paramagnetic) configuration. 
When the value of the exchange interaction Jh is large 
enough, crystal field effects, associated with the elonga- 
tion of the RuOe octaedra (A >0), tend to stabilize fer- 
romagnetic spin configurations. 

Indeed, for positive A, the dxy orbital is lower in en- 
ergy than the others two degenerate dxz.yz- 

Then, if the 

Coulomb repulsion is larger respect to A and the other 
parameters involved, one would occupy the first lower 
energy orbital and put the other hole in one of the two 
degenerate dxz,yz- In this case there is a gain in the ki- 
netic energy if the holes on the neighbour site have the 




0.50 



FIG. 4. Phase diagram for the two-sites cluster. Ferro, ferri 
and para stand for ferromagnetic, ferrimagnetic and param- 
agnetic ground state 



same spin of the hole moving. We have a kind of dy- 
namical double- exchange mechanism since the spin of the 
neighbour site is fluctuating in amplitude and phase, in- 
stead of having only phase fluctuation as happens in the 
manganite oxides. Still, strictly speaking, the exchange 
mechanism does not occur between localized spins. 

On the other hand, the gain in kinetic energy due to 
the large d — p hybridization, inverting the order of the 
molecular orbitals respect to the result of crystal field p^ ] 
and/or the compression of the RuOg octaedra, would pro- 
duce negative value for A and stabilize a paramagnetic 
configuration. 

Indeed, for negative A, the lowest energy configuration 
is realized by occupying with two holes the dxz,yz orbitals. 
In this case for large Coulomb repulsion we would have 
and effective Heisenberg antiferromagnetic Hamiltonian 
with a net exchange J ~ At'^/U which sets the lowest 
spin configuration in the ground state. 

These two limiting cases come in competition when | A | 
is small and of the order of t. In this case the gain in the 
kinetic energy allow to have configurations with non in- 
teger occupation number on each t2g orbital. This is the 
relevant situation which occurs in the ruthenium oxides 
either for the cubic and the tetragonal compound. Some 
comments on the limits of this cluster calculation are 
worthwhile. Within this two-sites problem one could not 
include any band effect which, for example, is relevant to 
understanding the dynamical nesting at incommensurate 
wave vector in the Sr2Ru04. Neverthless, the present cal- 
culation yields useful insights on the interplay between 
the magnetic correlations and the orbital degrees of free- 
dom in the complete series of the SrRuO compounds. We 
point out that the study of the electron dynamics in a 
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(P configuration have been attracted new interest also in 
other compounds hke the V2O3, suggesting more detailed 
analysis in this direction. 

In conclusion, in this paper we described the non- Fermi 
liquid behaviour observed in the normal state properties 
of Sr2Ru04. We have shown that, assuming that this 
compound lies very near to a quantum critical point, the 
temperature variation of various physical quantities is 
governed by spin fluctuations, whose effect is observed 
over a wide temperature range. A consistent explana- 
tion of the experimental results is obtained within this 
picture and assuming that only the electrons belonging 
to dxy band are involved in this mechanism. Moreover, 
we have given indication how to drive the system from 
a non-magnetic to a magnetic configuration both on the 
theoretical and experimental side which follows directly 
from the peculiar interplay between the lowering of the 
crystal symmetry, which involves the orbital degrees of 
freedom, and the spin correlations. 
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